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ABSTRACT 

Quasi-periodic oscillations and dwarf nova oscillations occur in dwarf novae and 
nova-like variables during outburst and occasionally during quiescence, and have ana- 
logues in high-mass X-ray binaries and black- hole candidates. The frequent low coher- 
ence of quasi-period oscillations and dwarf nova oscillations can make detection with 
standard time-series tools such as periodograms problematic. This paper develops tools 
to analyse quasi-periodic brightness oscillations. We review the use of time-frequency 
representations in the astronomical literature, and show that representations such as 
the Choi- Williams Distribution and Zhao-Atlas-Marks Representation, which are best 
suited to high signal-to-noise data, cannot be assumed a priori to be the best tech- 
niques for our data, which have a much higher noise level and lower coherence. This 
leads us to a detailed analysis of the time- frequency resolution and statistical properties 
of six time-frequency representations. We conclude that the wavelet scalogram, with 
the addition of wavelet ridges and maxima points, is the most effective time- frequency 
representation for analysing quasi-periodicities in low signal-to-noise data, as it has high 
time-frequency resolution, and is a minimum variance estimator. 

We use the wavelet ridges method to re-analyse archival data from VW Hyi, and 
find 62 new QPOs and 7 new long-period DNOs. Relative to previous analyses, our 
method substantially improves the detection rate for QPOs. 

Subject headings: methods: data analysis — techniques: photometric 



Introduction 



Cataclysmic variable (CV) light curves frequently contain multiple (q uasi-)periodic component s 
covering a wide frequency range, from 3 s dwarf no va oscillations (DNQs^ (IWoudt &: Warnerll2002al ) 
through ~100 s longer period DN Os QpDNQsl dWarner. Woudt fc Pretoriusl jooj) to ~5000 s 
quasi-periodic oscillations (QPOs) (jPatterson et alJll977l ). Each of these components can vary in 
period, amplitude and phase on short time-scales, against a high level of background noise. DNOs, 
with Q = ~ 10 4 — 10 6 , are generally more coherent than QPOs (Q ~ 1 — 10), but both are 
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seldom visible directly in the light curve, and the p eriodogram, which is desig ned for stationary 
phenomena, can fail to detect them (see for example iRobinson &: Natherlll979l ). Thus in order to 
study these and similar phenomena it is important to find a tool that detects them reliably. An 
investigation of time-frequency representations, which allow a non-stationary signal to be visualized 
in the time and frequency domains simultaneously, and have been successfully used in other areas 
of astronomy, would seem to be an intuitive and instructive place to begin. 

The most common time-frequency representation in use in astronomy is the wavele t trans- 
form, either in the form of the weighted wavelet Z-transform (WWZ) of iFosterl (|1996l ) or the 
wavelet scalogram. We will focus on the wavelet scalogram as our data are regularly spaced, 
and the WWZ is designed for irregularly spaced data and has poorer time-frequency resolution. 
The wavel et scalogram has been used to analyse, amongst other as tronomi cal phenomena, SR- 
type stars (iChinarqva &: Andronovl l2000 



Halevin et al 



2002 



Kiss. Szabo fc Beddinell200ffl CVs (Fritz fc Bruch 



2I) and planet a ry ne bulae ( Gonzalez Perez et al. 2003 ). The significance testin. 



1998 



method of iTorrence &; Compd (119981 ) for the wavelet scalogram is used mainly in solar astro- 
physics, where it h as been applied to magnetic fiel d measurements (jBoberg fc Lundstedtl I2002Q . 
coronal variability (IMarsh. Walsh Bromagd I2OO20 . high-fr equency oscillations in coronal loops 
(jKatsivannis et al.ll2002l ) and chromospheric UV oscillations (jMcAteer et al.l 12004 ) . 



Several other time-frequency representations have also been used in astronomy. The Gabor 



spect rogram has been applied to synthetic multi-co mponent signals and roAP stars (jBoyd et al 



995) and QPOs in b oth the atoll source 4U 1820-30 (jBelloni. Parolin fc Casellall2004h and Her X-l 



O'Brien et alJl200lh . The Choi- Williams distribution was applied to Wolf-Rayet stars (|Marchenko Moffat 



1998) and the Co ne-Kernel Representation to irregular pulsators ( Buchler. Kollath Cadmus! 
20041 ) and T UMi (jSzatmarv. Kiss fc Bebesll2003l ) . 



There are several key differences between our data and those analysed in the literature which 
make the choice of the most appropriate time-frequency representation for our analysis unclear. 
Firstly, in all of the papers mentioned above (and, indeed, all those in the literature) either visual 
inspection of the light curve or calculation of the periodogram is the primary source of periodicity 
detection; the time-frequency representation serves merely as an adjunct. For our signals, which 
have lower signal-to-noise ratios than those in the literature and which cannot be seen directly in 
the light curve or the periodogram, the time-frequency representation will be the primary means 
of detection, and hence it is imperative that we be able to test the significance of any detections 
statistically. Apart from iBovd et al.l (|1995l ). the only papers conduct i ng sig nificance testing of 
time-frequency representations are those using the ITorrence Compd (|1998l ) wavelet scalogram 
method. 

Our analysis also requires a broader frequency range than any discussed in the literature, as well 
as very high frequency resolution, particularly in the QPO frequency range as QPOs often appear 
as closely spaced harmonics. In addition, since both DNOs and QPOs are intermittent phenomena 
that can change their characteristics on a time-scale of 10 s of seconds, we need high time resolution 
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over the entire frequency range. Thus time-frequency representations with higher time-frequency 
resolution would be preferable for our analysis. Comparisons of the wavelet scalogram, Gabor spec- 
trogram, Choi- Williams distribution and the Cone-Kernel Representation show that, at least for rel- 
atively large amplitude variables stars such as irregular pulsators or Mira-types, the Choi- Williams 
distribution and Cone- Kernel Representation can give better time-frequency resolution than th e 



wavelet scalogram (e.g. iKollath k Buchlerlll997l : iBuchler k Kollath 



2001 



Kiss k Szatmarvl 12002 ). 



A plethora of papers investigate the effe ctiveness of different time- frequenc y representations us- 
i ng va ri ous synthetic signals (see for examplelGoupil. Auvergne k Baglin (119911) JCarothers k Swenson 



1992), 



1995), 



Sza tmarv. Vinko. k Ga] 



Buchler k Kollath 



( 1994 1) .iBovd et al 



(200 lj), iKiss k Szatmarvl ( 




Li. Jia k Zhao! (11 999) 



Hlawatsch et al 



Thavaparan k Kennedvl (|2003l ). ISzatmarv. Kiss k Beb 



20031 ) or iFigueiredo et al.1 (J2004J)), but as none of these synthetic data have features comparable 



to those of our data we cannot necessarily use their conclusions. 

Section [2] gives a brief introduction to the Wigner-Ville Distribution, and explains why the 
interference terms which distort this time-frequency representation have resulted in the develop- 
ment of a large number of alternative representations. In section [3] we analyse the time-frequency 
resolution and statistical properties of six of these alternative time-frequency representations using 
three synthetic (noise-free) signals mimicking key features of our data. We conclude that the Morlet 
wavelet scalogram, with the addition of wavelet ridges and maxima lines is the optimal choice for 
our data, both in terms of its time-frequency resolution and statistical properties. We use this 
representation to re-analyse archival data of VW Hyi in section [U and show that it does indeed 
provide a reliable tool for detecting QPOs, with a significantly higher detection rate than previous 
methods. Section [5] concludes. 



2. A brief introduction to time-frequency representations 

The archetypal time- frequency representation (TFR) of a signal x(t) is the Wigner-Ville Dis- 
tribution (WVD), given by 

/+oo 
x(t + T)x*(t-T)e 4 ^ Ti dT (1) 
-oo 



(|Villdll948l ). 



In practise, for finite duration signals, the pseudo Wigner-Ville distribution (PWVD), given 

by 

/+oo 
h(T)x(t + T)x*(t-T)e AnfTi dT. (2) 
-oo 

is used, where h{t) is a rectangular window that is zero outside the duration of the signal. If x(t) 
is sampled at times t = 0, At,2At, ... then, in order to prevent spectral aliasing in the WVD and 
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PWVD, either the signal must be oversampled by a factor of at least 2, or the analytic signal 
([Auger et al. Ill996l ). The analytic version x a {t) of a signal x(t) is the inverse Fourier transform of 



2x(f) if/>0 
if / < 



(3) 



We assume throughout the rest of the paper that the analytic signal is used. 

Ville originally investigated TFRs in order to study the instantaneous frequency of signals with 
time- varying components. A real (mono-component) signal x(t) may be written as a time varying 
amplitude a(t) modulated by a time varying phase 4>: 



x(t) = a(t) cos[(p(t)} 



(4) 



with a(t) > 0. The instantaneous (angular) frequency is then defined as the positive derivative of 
the phase: 

u(t) = 4>'{t) > (5) 



(jGoncalves et al.lll997l ). Since there are many possible choices of a(t) and <p{t) for any given signal, 
uj{t) is not uniquely defined. However, using the analytic signal, we have x a (t) = a(t) exp[ic/>(i)] and, 
since x(t) = Re[x a (t)}, then x(t) = a(t) cos 4>{ t). We call a(t) the analytic amplitude of x(t) and 
4>'{t) its instantaneous frequency (|Mallatlll998l ). Multicomponent signals in which the components 
occur in disjoint frequency bands can be thought of as the sum of mono-component signals; it is 
the aim of the TFR to show the instantaneous frequency of each component accurately. 

Figure QJa) shows the time-series of a synthetic signal which mimics a noise-free dwarf nova 
(DN) lightcurve sampled with an integration time of 4 s, with a first harmonic DNO changing 
linearly in period from 25 s to 30 s over the duration of the run, fundamental and first harmonic 
QPOs starting with periods 750 s and 375 s respectively and tracking the period change of the 
DNO, and an lpDNQ wi th a constant 70 s period, based on run s0484 of VW Hyi, discussed in 
Woudt fc Warnerl (|2002al ). An idealised TFR (Figure C^b)) and the PWVD (Figure Q^e)) are shown 
below the time-series. 

While the changing DNO and QPO periodicities are clearly visible in the PWVD in FigureQTc), 
as is the constant period lpDNO, there are also five sets of spurious signals, called interference terms 
(ITs), which are an artefact of the quadratic nature of the WVD. The WVD of an N-componen t 



signal will always consist of N signal terms and N ^ ^ ITs dMecklenbrauker fe Hlawatschl Il997l ) 



ITs oscillate in the time and frequency directions (jHlawatsch &: Urbankd Il994l ) , and the closer 
together two components are in freque ncy, the slower their ITs will oscillate in the time direction 
(jHlawatsch &: Boudreaux-Bartels!ll992l ). This is clear in figure (H|c) : the ITs due to the interaction 
between the two close QPO terms oscillate much more slowly in the time direction than the ITs 
due to the other interactions, whose generating signal terms are further apart in frequency. 

The presence of ITs in the WVD make it unsuitable for analysis of data in which the true 
signal components are unknown, or cannot be inferred by other means - as is the case for our data. 
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Fig. 1. — (a) Time-series of a synthetic signal which mimics a noise- free dwarf nova (DN) lightcurve, 
with a first harmonic DNO changing linearly in period from 25 s to 30 s, fundamental and first 
harmonic QPOs starting with periods 750 s and 375 s respectively and tracking the period change 
of the DNO, and an lpDNO with a constant 70 s period, (b) Idealised TFR of the signal in Figure 
1(a). (c) WVD of the signal in Figure 1(a) 
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However, since the ITs are oscillatory, they can be reduced or removed entirely by smoothing the 
WVD in the time and frequency dir ections. IT attenuation is usually ach ieved by convolving the 
WVD with a 2-D smoothing kernel (jHlawatsch &: Boudreaux-Bartelslll992l ). 



IT reduction comes at a cost: t he more smo othing in time and/or frequency, the poorer the 
resolution in time and/or frequency (ICohenl ll989). The WVD has the best possible time- frequency 
resolution of all TFRs: for a time-dirac signal centered at to, W x (t, f) = 5(t — to) , and for a 
frequency dirac signal centered at /o, W x (t, f) = 5(f — fo) (|MaHatlll998l ). All other TFRs generally 
have poorer time and frequency resolution, or may achieve that of the WVD for certain types of 
signal only. 

Since the geometry of the ITs depends on the components of the signal, different smoothing 
kernels have been developed for different kinds of signals: there is no single kern el which is optimal 
for all signals. Most ke rnels however belong to the Cohen class (|Cohenlll989l ) or the affine class 
(jRioul &: Flandrinlll992l ) or both, and have one or more parameters which can be varied to choose 
the amount of smoothing in time and/or frequency. 



The Cohen TFRs (ICohenlll989l ). which we will denote generally as C x (t,f), are derived from 
a time- frequency convolution of the WVD with a smoothing kernel ip(t, /): 



C x (t, f)= t', f - f')W x (t', f')dt'df 

J t > Jy 

or, more traditionally, as 

C x (tJ)= [ [ *(T,v)A x (r,vy 2 ^-^dTdv 

J T J V 



(6) 



(7) 



where ^(r, v) is the Fourier transform of ip{t, f) and A x (t, v) is the ambiguity function (AF) 
of x(t), defined as: 



K x {t,v) =2j t x(t + T)x*(t-r)e 



dt 



2j f x{f + u)x*{f-u)e l4 ^fdf. 



(8) 



Here r is the time-delay, and v is the Dop pler value, which meas ures the difference in fre- 
quencies between points at a given time-delay (jHlawatsch et al.lll995l ). The AF, which is the 
Fourier transform of the WVD, can b e interpreted as a joint time-frequency correlation function 
(jHlawatsch Boudreaux-Bartelslll992l ). 



Each member of Cohen's class is thus associated with a unique, signal independent kernel 
function ^(r, u); the WVD is a member of the Cohen class, with kernel ^wvd( t , v) = 1. 



The general affine TFR (A x (t,a;TL)) of an analytic signal x(t) is given by 
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A x (t,a;U)= / n( ,at)W x (8,Z)dsde (9) 

«/ — oo J —oo " 

where z/) is a smoothing kernel. Equivalently, using the AF instead of the WVD, we have 

/OO /"OO 
/ ^{a^T/a^d^e-^didT (10) 
-oo J —oo 

where ^(^, r) is the ambiguity domain version of II(£, v). Affine TFRs are time-scale represen- 
tations rather than time- frequency representations, but a acts as a measure of th e frequency: reduc 



inga reduces the time support of Ht, s , and hence increases its frequency range ([Goswami fc Chan 



19991 ) . Since the WVD is scale-invariant, it is also a member of the affine class. 



In both Cohen and affine TFRs, the kernel determines the amount of smoothing in time and/or 
frequency. The important difference between affine and Cohen TFRs is that the affine kernel is 
scaled depending on the analysing frequency, while Cohen kernels retain the same 'size' at all 
frequencies. This means that for affine TFRs the time and frequency resolutions depend on the 
frequency at which the signal is being analysed, while for Cohen TFRs the time and frequency 
resolutions are the same at all frequencies. 



3. Six smoothing kernels investigated 

We will investigate 4 Cohen and 2 affine smoothing kernels which have either been previously 
used in used in astronomy or have properties which suggest they are appropriate for our data; Table 
[T] gives a summary of the main features of the kernels discussed. We begin our analysis by finding 
the optimal parameters needed for each TFR in order to represent the synthetic DN lightcurve 
introduced in section [2] without interference terms. This signal tests the frequency resolution of 
each TFR over a broad frequency range, and ensures that the values of the parameters chosen 
completely attenuate interference terms over the whole frequency range. Using these optimal 
parameters, we then further investigate the time resolution of each TFR by using intermittent 
sinusoids of different frequencies. Where possible, we include a discussion the use of the kernel in 
the astronomical literature. We conclude this section with a discussion of the statistical properties 
of the six kernels. 

3.1. The Smoothed Pseudo WVD 

One method of smoothing the WVD is to use the separable kernel 



T T * 

^spwvd(t, v) = g(- )g*(- -)h(-u) 



(11) 
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Table 1. TFR kernels and their properties 



Name 



Kernel \I/(t, v) Parameter (s) Cohen Affinc 



Statistics 



Wigner-Ville (WVD) 



1 



Snoise{f^) 2 



Smoothed Pseudo g{^)g*{-^)h(v) 

Wigner-Ville (SP- 

WVD) 



Lg, Lift 



2M—1 '-'noise 



(/)¥ 



Choi- Williams (CWD) 
Born- Jordan (BJD) 
Cone-Kernel (CKR) 



sin(7rri^) 

7TTIS 



/ 



Snoise (/) 2 2 



Snoise (/) 2 2 



^noise \J ) 9 



Affinc Smoothed 
Pseudo Wigner-Ville 
(ASPWVD) 



2.u. 



2 _<? - f^Xa. 

1 '-'rocnse V* 3 J 2 



Gabor Spectrogram Wh{—T,—u) 
(SPEC) 



a . ( f\X2 

u noise \J ) 2 



Scalogram (SCALO) 



$noise{f) 2 
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with time-smoothing window g(t) and freq uency-smoothing win dow h{v) (which is the Fourier 



transform of a time-smoothing window h(t)) (jAuger et al. II1996I ). The length L g (L^) of the 



window g(t) (h(t)) controls the amount of smoothing in the time (frequency) domain. A longer 
g(t) gives more time smoothing and hence worse time res olution, and a longer h(t) gives less 



frequency smoothing and hence better frequency resolution (jHlawatsch et al.lll995l ) 



Figures [2] and [3] show the smoothed pseudo WVD (SPWVD) of the synthetic DN lightcurve, 
with L g = 200 s in both figures, but Lh = 200 s in the first, and Lh = 1200 s in the second. 
The effect of increasing Lh (i.e. decreasing the frequency smoothing) is marked: in the first figure, 
we cannot differentiate the two long period QPO components, but in the second they are clearly 
visible. Interference terms between the two long-period components in Figure [3] are still prominent. 
By increasing L g to 1000 s, as in figure^ the interference terms are completely attenuated. Notice, 
however, the effect of the extensive time smoothing at the ends of the DNO. 

Next, we use these optimal values of L g and Lh to investigate an intermittent sinusoid of period 
100 s, shown figure E^a), mimicking for example an intermittent lpDNO. The integration time is 
4 s. 

The time-smoothing introduced by increasing L g , which was necessary to attenuate the ITs, 
has come at a price: the ability of the SPWVD to clearly resolve temporal discontinuities in the 
signal. In figure EJ^b) the smoothing in the time domain results in a SPWVD that does not clearly 
show the beginning or end of each sinusoidal burst. Indeed, it gives the impression that the sinusoid 
is perhaps present for the duration of the signal, but at varying amplitude. There are two reasons 
for this: signal terms, which should not be smoothed, are smeared in the time dimension by the 
long time window, and ITs which oscillate in the frequency direction between the signal terms are 
not removed because we have set the length of the frequency smoothing window to ensure that we 
have optimal frequency resolution (i.e. minimal smoothing). 

We are thus presented with the central issue that dogs most Cohen kernels: if ITs are attenuated 
by smoothing in time, the result is poor time resolution, which is particularly problematic at high 
frequencies. 



3.2. The Affine Smoothed Pseudo WVD 

The affine smoothed pseudo WVD (ASP WVD) is the affine counterpart of the SPWVD, and 
its kernel is given by 

Vaspwvd = h(-)g{- — -), (12) 
a a 

where g(t) and h{t) are smoothing windows as for the SPWVD. The lengths of the time and 
frequency smoothing windows, L g and Lh respectively, again allow a flexible choice of time and 
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Fig. 2. — Smoothed Pseudo Wigner-Ville distribution of synthetic DN lightcurve, with L g 
L h =200. 
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Fig. 3. — Smoothed Pseudo Wigner-Ville distribution of the synthetic DN lightcurve, with 
L g =200s, L h =1200s. 
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Fig. 4. — Smoothed Pseudo Wigner-Ville distribution of the synthetic DN lightcurve, with 
L q =1000s, L h =1200s. 
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Fig. 5. — Smoothed Pseudo Wigner-Ville distribution of an intermittent sinusoid of period 100 s, 
with L g 1000s, L h 1200s. 
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scale (frequency) resolution. 

Figures [6] and [7| show the ASPWVD of the synthetic DN lightcurve, with L g set to s in the 
first case (no time smoothing) and 320 s (at the largest scale) in the second. is 320 s in both 
cases, which is the minimum length needed to resolve both QPOs. What is apparent in Figure 
[7J when compared with the SPWVD shown in Figure HI is that less time smoothing is needed to 
remove the ITs in the ASPWVD, resulting in less noticeable end-effects. However, the frequency 
resolution decreases in the ASPWVD as the frequency increases, so the DNO in the ASPWVD is 
broader than the DNO in the PSWVD. 

In comparison with the SPWVD, the ASPWVD of the intermittent sinusoid, shown in figure[8j 
has good time resolution (since for affine distributions the time resolution improves as the frequency 
increases), but has poorer frequency resolution than the SPWVD. No ITs are visible between the 
bursts of signal. 

The top panel of figure [9] shows a signal mimicking an intermittent 20 s DNO with changing 
amplitude, sampled with an integration time of 4 s. The ASPWVD of this signal, calculated 
using L g = Lh = 320 s is shown in the bottom panel. Very broadly smoothed signal terms and 
ITs considerably distort the true signal, and it is not possible to remove these interference terms 
without losing the optimal resolution achieved for the analysis of the DN signal. 



3.3. The Choi- Williams Distribution 



Instead of using a separable smoothing kernel, several kernels have been developed that depend 
on the product of r and v. One of the earliest examples is the Choi- Williams distribution (CWD) 
(|Choi &: Williamslll989l ). which is characterized by the kernel 



CWD 



{r,y) 



(2tttv) z 



(13) 



a is a positive paramet er that jointly controls smoothing in time and frequency: a larger value 
for a gives less smoothing (jHlawatsch et al.lll995l ). In fact, as a — > +oo, the CWD tends to the 
WVD. Pure sinusoids have the ideal time- freq uency resolution of the WVD , but all other signals 
are broadened when compared with the WVD (jHlawatsch fc Urbankdll994l ). 



Figure [10] shows the CWD of the synthetic DN lightcurve with a = 100. Reducing a to 0.1, 
as shown in figure [TT| which increases the amount of smoothing, does not fully eliminate the ITs; 
further reductions in a have no noticeable affect. In fact, ITs due to different frequ ency components 



occur ring at the same time c annot be attenuated completely if we use the CWD (jHlawatsch et al 



19951 ) . In iKiss &; Szatmaryl (j2002j) , the good time resolution of the CWD enables the increase in 
frequency of the main pulsation to be seen clearly, and the strong interference terms caused by the 
interaction of the main component with its harmonics are not a problem, since they can be easily 
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Fig. 6.— ASPWVD of the synthetic DN lightcurve, with L g = 0s and L h = 320s. 
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Fig. 7.— ASPWVD of the synthetic DN lightcurve, with L g = 320s and L h = 320s. 
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Fig. 8. — ASPWVD of an intermittent sinusoid of period 100 s, with L„ = 320s and Lh = 320s. 
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identified as spurious using the periodogram. However, for our data, since complete IT attenuation 
is required, the CWD is not ideal. 



3.4. The Born-Jordan Distribution and Cone-shaped Kernel Representation 

If we require a TFR that preserves finite time- and frequency-support, then the simplest choice 
for ^(r, v) is 

sin(7rrz^) 

^bj{t,u) = — -, (14) 

TTTV 

which defines the Born-Jordan distribution (BJD). This kernel does not admit smoothing of ITs 
in the frequenc y direction. Allow ing smoothing of the BJD in frequency gives the Zhao- Atlas-Marks 



representation (jZhao et al.lll990l ). also known as the cone-shaped kernel representation (CKR), with 
the kernel 

, . . ,sin(7rr^) . . 

®CKr(t,v) =g( T )\T\—± (15) 

TTTV 



where g(r) is any window function (jHlawatsch et al.lll995l ). The length of the window, L 



again controls the amount of frequency smoothing. 



Figure [T2l shows the CKR of the synthetic DN signal with L g = 1000 s, which is the minimal 
length required to completely attenuate the ITs between the QPOs. The frequency resolution is 
comparable to the SPWVD and ASPWVD, and the edge effects are less marked than those of the 
ASPWVD. Notice, however, that the way in which the CKR smooths the signal in figure [T2l causes 
the QPOs to appear to oscillate slightly in frequency. IT attenuation in the CKR is complex. 
For components which are close in frequency, the attenuation will depend on the time distance, r, 
between the components. If the time distance between two components is r, and < |r| < \L g 
then the resulting IT will be amplified, rather than attenuated. However all terms (signal and IT) 
with \t\ > ^L g will be completely suppressed. We are thus faced with a dilemma: increasing L g 
improves frequency resolution, but also makes it more likely that ITs will fall in the < |r| < \h g 



range, and hence be enhanced (jHlawatsch et al.lll995l ). 



The CKR of the intermittent 100 s is shown in figure [TBI In comparison with the previous two 
TFRs, the frequency resolution of the CKR is excellent. The time resolution is similar to that of 
the ASPWVD. ITs, although of low amplitude, are present, again appearing as a low amplitude 
100 s oscillation between the true signal bursts. 

Figure Q3] shows the CKR of the intermittent, amplitude-modulated 20 s oscillation signal. 
The amount of smoothing in time required to remove the low-frequency ITs in the synthetic DN 
lightcurve results in the signal terms of the intermittent DNO being completely smoothed. In 



Fig. 9. — ASPWVD of an intermittent, amplitude modulated sinusoid of period 20 s, with L g = 320s 
and L h = 320s. 




Fig. 10 — 



CWD of the synthetic DN lightcurve with a = 100. 
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addition, the ITs are amplified (due to the complicated IT attenuation of the CKR kernel) and 
smoothed , resulting in the appearance of five continuous frequency bands, even more pronounced 
than those in the ASPWVD of the same signal. While this is a somewhat extreme example, it does 
show that the effects of the CKR must be carefully monitored. 

As the CKR can attenuate signals in a manner that is not easy to predict, and can amplify 
ITs rather than attenuating them, it is not ideal for analyses requiring significance testing. For 
large amplitude signals where this is not the case, the CKR can be an extremely useful tool: in 



figures 2 to 4 of iBuchler. Kollath Cadmusl (|2004l ). for example, the CKR is seen to give sharp 



images for the five large-amplitude irregular pulsators under discussion. The fundamental period 
and harmonics are visible; and since these are relatively widely spaced in frequency, the ITs are 
easily identifiable. 



3.5. Gabor Spectrogram 

The (Gabor) spectrogram's kernel is given by 



^SPEC = W h {-T,-v), 



(16) 



where Wh(—T, —u) is the WVD of h(t), any short-time analysis window localised around t = 
and / = 0. This means that the time-frequency resolution of the spectrogram is dictated by the 
resolution of h(t), rather than that of the analysed signal. Smoothing in the spectrogram is thus 
extensive, r esulting in virtually complete interferen ce attenuation, but also in poor time-frequency 



resolution (jHlawatsch 8z Boudreaux-Barteld Il992l ). The choice of h(t) determines the trade-off 



between the time spread and the frequency spread of the smoothing function: time resolution 
is proportional to the time duration of h(t), while frequency resolution is proportional to the 
bandwidth of h(t). 

Figure [T5l shows the Gabor spectrogram of our synthetic DN signal, using a Gaussian window of 
length 2000 s, which is the shortest that allows resolution of the QPO components. Each component 
is clearly visible, and there are no interference terms, although some smoothing is apparent near 
the ends. 

Figure [16] shows the Gabor spectrogram of an intermittent sinusoid, using the same window 
as in the previous example. We se e that requiri ng go od frequency reso lution has again resulted 
in poor time resolution. Indeed, as ICohenl (j 19891 ) and IZhao et al.l ()1990l ) point out, if one wishes 
to have accurate time and frequency measurements using the spectrogram, two separate analyses 
must be undertaken, using windows of different length. 



_ Boyd et al.1 ()1995l ) use the Gabor transform to investigate HD 60435 (see fig. 7 of lBoyd et al . 
Because they are analysing a specific oscillation (11.64 min) they are able to choose an 
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Fig. 11. — CWD of the synthetic DN lightcurve with a = 0.1. 




Fig. 12.— CKR of the synthetic DN lightcurve with L g = 1000s. 
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Fig. 13. — CKR of a intermittent sinusoid of period 100 s, L 9 =1000s. 
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analys is window width which gives the required time-frequency resolution in this narrow frequency 
band. iBelloni. Parolin <fe Casellal (|2004l ) also consider a relatively narrow frequency b a nd, a nd are 
hence able to achieve an appropriate time-frequency resolution balance. iDolan et al.l (|1998l ) gives 
one of the the few examples of the analysis of a transient signal with the Gabor spectrogram, but 
since the frequency range is relatively narrow, they are able to choose a window length which gives 
the optimal time-frequency resolution for analyzing the signal. 



3.6. The Wavelet Scalogram 



The wavelet scalogram is the affine counterpart of the Gabor spectrogram. Just as the Gabor 
spectrogram can be calculated using any short-time analysis window that is localised around t = 
and / = 0, the wavelet scalogram can be calculated using any "mother wavelet" V'o that is centered 
in the neighbourhood of t = 0, has zero av erage (J°° ipo(t)dt = 0) and is normalized to have unit 
energy ( \ipo(u))\ 2 duj = 1) ( Mallatl 1 19981 ). From the mother wavelet, ipo, we define a family of 
functions 



Vv( s ) 



i 



t 



s > 0, b e 



(17) 



Then for a given ipt,a, the wavelet scalogram, Ascalo, has the affine kernel 



^scalo{t,v) = Ww, t Jr,u 



(18) 



where W^ t o (r, v) is the WVD of ip tja . 

The choice of mother wavelet depends on the signal to be analysed and the type of analysis 
desired. If the shape of the signa l to be detected is known, the n the wavelet with the nearest 
matching shape should be chosen (jMevers. Kelly fc Q'Brienlll993l ). Considering the often (ampli- 
tude modulated) sinusoidal shape of DNOs and QPOs, smoothly oscillating wavelets such as the 
Morlet, Mexican Hat and Paul wavelet are indicated, rather than the abrupt step function of the 
Haar wavelet, or the spiky Daubechies wavelet. 

An additional consideration is that in the scalogram of a complex wavelet, such as the Morlet, 
phase information from the data is lost. A scalogram using a complex wavelet has constant power 
across the time duration of an oscillation, and the time locations of maxima and minima in the 
oscillation cannot be detected. For real wavelets, however, such as the Mexican Hat wavelet, 
phase and amplitude information i s superposed in the time-frequency plane and extrema are easily 
detected (jPercival fc Waldenl |2000j) - a feature that is not present in any of the TFRs studied so 
far in this paper. 



In order to explore fully the behaviour of the wavelet scalogram we thus investigate two smooth 
wavelets with very different analytical capabilities: the complex Morlet wavelet and the real Mexi- 
can Hat wavelet. 
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Fig. 14. — CKR of an intermittent, amplitude modulated sinusoid of period 20 s, with L g = 1000s 
and L h = 1000s. 
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Fig. 15. — Gabor spectrogram of the synthetic DN lightcurve, with L g = 2000s. 
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Fig. 16. — Gabor spectrogram of an intermittent sinusoid, with window length 2000s. 
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3.6.1. The Morlet Wavelet 



The Morlet wavelet is described by 



7T 



I _ I*,2 

4 e 2 



(19) 



where ?? = i/s is a non-dimensional time parameter and wq is the wavenumber (jTorrence fc Compo 



19981 ). The Morlet wavelet is ess entially a complex expon ential of frequency modulated in am- 



plitude by a Gaussian window (jPercival k, Waldenl 12000 ) . We have chosen to use wo = 6 in our 
analysis, which gives a Morlet wavelet which has 5 oscillations. This is the most commonly used 
value of cjo, but is also appropriate given the QPO detection criteria which will be discussed in 
section HI 

Since the wavelet is localized in time and scale, the wavelet scalogram at a particular time t' 
is affected only by signal values close to t' , in a radius that depends on the width of the wavelet 
as determined by the scale a. The cone of influence (COI) of a point (t',s) on the time-scale 
plane is the set of all points (t",s) for which x(t') affects Ascalo{^, a). Finite data length means 
that for points near the beginning or end of the data set, the support of AscALo(t, o.) effectively 
includes a number of points with zero amplitu de (i.e. points before or aft er measurement started) , 
and therefore has significantly reduced power. iTorrence &: Compd (|1998l ) define the e- folding time 
t s for the wavelet as the time-radius at which the wavelet power first drops by a factor of e~ 2 
at an edge. The time extent to which these 'edge effects' affect the wavelet scalogram depends 
on the scale of the wavelet and the type of wavelet; at shorter scales (shorter periods) the effect 
is minimal, while at very long scales, the effect is marked. For the Morlet wavelet, t s = \/2s 
(jTorrence Compolll998l ). 



It is preferable to plot the wavelet scalogram as a function of frequency rather than scale. For 
the Morlet wavelet, the equivalent Fourier period for a given scale s is 



Period 



4ns 



u 



+ V2 + ojq 1 



(20) 



([Meyers. Kelly Q'Brienlll993l ; ITorrence fc Compolll998l ). Figure [T71 shows the Morlet scalo- 
gram of the synthetic DN lightcurve. While all components are clearly visible, and there are no 
interference terms, the lpDNO and especially the DNO are not well resolved in frequency. The large 
black 'U' marks the e-folding time at each scale: points outside this region are distorted by edge 
effects. It is clear that while edge effects are minimal at high frequencies, they become significant 
at low frequencies. 

In both the Morlet scalogram of the intermittent sinusoid with period 100 s (figure [T8|) and the 
intermittent, amplitude-modulated sinusoid with period 20 s (figure [T9]) . we have an IT-free TFR 
in which the signal is well resolved in time, but this has come at the expense of poor frequency 
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resolution: the Morlet scalogram has by far the worst frequency resolution of any of the TFRs 
investigated so far. There is, however, a way we can compensate for this poor frequency resolution, 
while maintaining the excellent time resolution and IT-free properties of this TFR. 



Delprat et al.l (|1992l ) show that for complex wavelets such as the Morlet, the wavelet scalo- 



gram of a monocomponent signal with slowl y chang i ng low frequency reaches a maximum at the 
instantaneous frequency (f>'(t) of the signal. iMallatl (|1998l ) extends this proof to include the ac- 
curate measurement of the instantaneous frequency of multiple components changing rapidly at 
high frequency. Thus the points (i,f) at which the Morlet scalogram AMORLEr(t, /) has a local 
maximum in frequency, i.e. for which 



9A M QRLET{tJ) 

df 



and 



d 2 A M QRLET{tJ) 



<o, 



(21) 



c an be used to plot the instantaneous frequency of our signal accurately in the form of wavelet 
ridaes dMallatl 1 1 9981 ). 



Figures [2"0j 1211 and 1221 again show the Morlet scalograms of our three test signals, but with the 
inclusion of wavelet ridges. Figure [T7] looks pleasingly close to the ideal TFR shown in figure [TJ 
Both figures l2"T1 and 1221 show wavelet ridges at frequencies other than the instantaneous frequency of 
the signal, but as these occur only at very low amplitudes they do not present the interpretational 
challenge that interference terms did in previous TFRs. 



3.6.2. The Mexican Hat Wavelet 



One of the most commonly used real wavelets is the Mexican Hat Wavelet (MHAT), given by 
the normalized second derivative of a Gaussian function: 
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Fig. 17. — Morlet scalogram of the synthetic DN lightcurve. 
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Fig. 18. — Morlet scalogram of the 100 s intermittent sinusoid. 
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Fig. 19. — Morlet scalogram of the 20 s intermittent, amplitude modulated sinusoid. 




Fig. 20. — Morlet scalogram of the synthetic DN lightcurve, with ridges. 
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Fig. 21. — Morlet scalogram of the 100 s intermittent sinusoid, with ridges. 
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Fig. 22. — Morlet scalogram of the 20 s intermittent, amplitude modulated sinusoid, with ridges. 
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^(r?) = (1 - 4r)e^ (22) 

where again rj = t/s. The equivalent Fourier period for a given scale s is 

Period = -JL. (23) 
a/2T5 



which effectively increases the minimum observed time-scale from 2c£i to 7.95dt (jMcAteer et al 



2007). As a result, the MHAT wavelet is not useful for examining very short period phenomena. 
Figure [23] shows the MHAT scalogram of the synthetic DN signal. Notice that the highest frequency 
achievable is 0.03 Hz, which is below the starting frequency of the DNO. It is clear that the frequency 
resolution of the MHAT wavelet is considerably worse even than that of the Morlet wavelet. The 
usefulness of the MHAT wavelet is, however, apparent in the appearance of the lpDNO and the 
QPOs in the scalogram: individual extrema are discernable (the wavelet scalogr am is an ent irely 



positive TFR, so all signal extrema appear as maxima in the wavelet scalogram (|Mallatlll998l )). 



Similar to the wavelet ridges of the Morlet wavelet, the points (t, f) at which the MHAT 
scalogram AMHArifi /) ^ as a l° ca l maximum in time, i.e. for which 

dA MH Mt,f) = Q and &A M HA T {tJ) < Q) (24) 



form maxima lines which identify maxima across frequencies ( MallatJ 19981). Maxima line s 



traditionally u s ed to detect singularities and discontinuities. iPolvgiannakis et al. I (|2003l ) and 



arc 



McAteer et al.l (|2007l ). for example, use maxima lines to investigate the sunspot index and so- 
lar flare x-ray emission respectively. For our data, however, maximal lines can be used to detect 
the extrema in quasi-periodicities, which can be used to measure coherence, investigate phase re- 
lationships between different frequency phenomena and construct O-C diagrams. Figures [Ml and 
[25] show the MHAT scalogram of the synthetic DN light curve and the 100 s intermittent sinusoid 
respectively, with the addition of maxima lines. 

Both the Morlet scalogram with the addition of wavelet ridges and the MHAT scalogram with 
the addition of maxima lines provide valuable information in the analysis of our data. To facilitate 
easy analysis, we combine the information provided by the wavelet ridges and maxima lines on a 
single 'Morlet ridged scalogram': we plot all wavelet ridges, but only those points of the maxima 
lines which intersect with wavelet ridges. Figure [26] shows the our final TFR of the synthetic DN 
signal. Here we finally have a TFR in which all frequency components can be seen clearly, as well 
as the times at which each component has an extremum. 

At least in the case of noise- free signals, it is clear that the Morlet ridged scalogram is the 
TFR best suited to our needs, providing excellent frequency resolution over a broad frequency 
ranges, combined with the time resolution needed to resolve intermittent signals. However, when 
analysing real data the statistical properties of the TFR are important for judging the likelihood 



that a detected oscillation is a true signal, and not an artefact of the noise. It thus behoves us to 
investigate the statistical properties of each of our six kernels, before declaring a winner. 



3.7. Statistical Properties of TFRs 



We assume that we have a non-stationary deterministic signal superimposed on weakly sta- 
tionarjQ non-white noise. This is not an unreasonable assumption: the noise spectrum of our data 
changes very slowly if at all over the duration of an observing run. We will discuss methods of 
estimating the power spectrum of the noise in section 14.2} but note here that the assumption of 
weak station arity ensures that the per iodogram can be used as an estimator for the power spectrum 
of the noise (jTimmer &; Koenielll995l ). 



The WVD of a re alization of a stochastic signal is an estimator for the theoretical WVD of 
the generating process (jMecklenbrauker fc Hlawatschlll997l ). For a weakly stationary noise process 
with power spectrum S(f), the expected value of the WVD, E[W x (t, /)], is given by the power 
spectrum of the process: 

E[W x (t,f)] = S(f) (25) 



(|Qiulll993l ). For a deterministic signal corrupted with independent, identically distributed (i.i.d) 
gaussian noise with power spectrum S no i se (f), the expected value of the WVD is the sum of the 
WVD of the signal and the spectrum of the noise: 



E[W x (t,f)]=W signal + S noise (f) 



(26) 



(jFlandrin &; Martini 1 19841 ). The WVD is thus a biased estimator of the signal, and the larger 
Snoise(f), the worse the bias. The WVD of the noisy signal, since it is determined by the power 
spectrum of the noise, is also not consistent: W x (t, /) ~ S no i se (f)-^-. 

Thus the WVD suffers from the same estimation problems as the periodogram. However, 
by estimating the power spectrum of the noise, confidence contours (the 2-D equivalent of confi- 
dence levels) can be used to determine regions in the time - frequency plane at which the signal has 
power significantly above that expected for pure noise (see iTorrence &: Compd (|1998l ) for a detailed 
discussion of confidence contours for the wavelet scalogram). 

The expected values of any Cohen TFR of a deterministic signal corrupted with i.i.d gaussian 
noise are also given by the sum of the noise-free TFR and the power spectrum of the noise, as was 
the case for the WVD. Thus these TFR estimators, which include the CWD, the BJD and the 
Gabor spectrogram, are also biased, with the amount of the bias depending on the noise level. We 



A time series is weakly stationary if all realizations {x t } of the process {Xt} have the same mean, the 
(auto)covariance between any 2 samples from the same reali zation is a function of the time-lag between the samples 
only, and the variance is finite (|Kendall, Stuart fc Ordll983h . 
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Fig. 23. — MHAT scalogram of the synthetic DN lightcurve. 
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also have C x (t, f) ~ S„ 



in a lower variance ( 


Williams 


(Hcaron & Amin 


1995), while 



Williamsl (|1998l ) shows that the Gabor spectrogram, when a rectangular window is used, gives the 
best performance of all Cohen TFRs in the presence of non-white noise, as it has the smallest 
variance. 

For the SPWVD with L„ = M, for a deterministic signal corrupted by i.i.d gaussian noise, 



SnoiseifW dFlandrin fc Martini 11984 1. Thus the variance of the SPWVD 



CsPWVD(t, f) ~ 2M-l '-' nmse ^ 

decreases as the length of the time-smoothing window increases. 

For the ASPWVD, with L g at scale s given by M s , then for a deterministic signal corrupted 
by i.i.d gaussian noise, 



AASPWVD(t, s) 



2M„ 



1 



(27) 



(jFlandrin Martinlll984l ) . Thus the variance of the ASPWVD decreases as the window length 
increases, as in the case of the SPWVD, but also depends on the analysing scale s. 

Since the CKR's integrals with respect to frequency, and over the entire time-frequency plane, 
are zero (because ^ckr( t , v) is zero on the entire z^-axis), the CKR is not an energy distribution 
i n the statistical sense (which is why we refer to it as a representation, rather than a distribution) 
(jHlawatsch et al.lll995l ). This means that its expected value is not related to the power spectrum 
as for the other TFRs. The CKR of a deterministic signal corrupted by i.i.d gaussian noise is the 
sum of the CKR of the noise-free signal and D(t, /; g), where 



D(t,f;g) 



g(T)\T\r(T)e- i27TfT dT 



(28) 



and r(r) is the auto-covariance of th e noise. This mean s that for wh ite noise, D(t, f; q ) = 



and the CKR is an unbiased estimator (|Qh &: Marks! Il992l ). However, iHearon &: Amid (| 19951 ) 
show that of all the Cohen TFRs the CKR has the highest variance, and hence the least desirable 
performance in the presence of noise. In addition, confidence contours cannot easily be constructed 
for the CKR since the stochastic distribution of the CKR is not calculable from the power spectrum 
as is the case for the other TFRs. 



Percivall (|1995l ) shows that the time integral of the wavelet scalogram can be used as an 
unbiased, consistent estimator for the power spectrum. It is also the minimum affine variance 
estimator, and we have Ascalo ~ \S n oise)di- Thus from a statistical point of view, the wavelet 
scalogram is also the best TFR to use for our data. 
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Fig. 24. — MHAT scalogram of the synthetic DN lightcurve with ridges. 
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Fig. 25. — MHAT scalogram of the 100s intermittent sinusoid with ridges. 
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Fig. 26. — Morlet ridged scalogram of the synthetic DN lightcurve with wavelet ridges and maxima 
points. 
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4. Application of Time-Frequency Techniques to VW Hyi 

In this section we use the Morlet ridged scalogram introduced in the previous section to analyse 
quiescent and outburst lightcurves from VW Hyi. Our purpose is threefold. Firstly since all of 
these runs have been analysed in detail previously we see how previously detected DNOs and QPOs 
appear in the wavelet scalogram, and hence build up a set of criteria for objectively detecting QPOs 
using the wavelet scalogram. We then use these criteria to see if there are QPOs in the VW Hyi 
data which have not been identified by the old methods. Finally we use the wavelet methods to 
investigate the non-stationary behaviour of all DNOs and QPOs detected. 

This analysis establishes the reliability of the Morlet ridged scalogram (we find that all previ- 
ously identified DNOs and QPOs are detected by the scalogram), and provides an opportunity to 
compare the information elicited by the new techniques with that given by traditional techniques 
such as the Fourier Transform and O-C curves. 



4.1. Observations 



VW Hyd ri is an SU UMa dwarf nova with an outburst time scale of the order of 28 days 
(|Warneiill995l ). Since it is circumpolar at the Sutherland site of the South African Astronomical 
Observatory, and is relatively bright (9-13 mag), VW Hyi presents an extended season for maximum 
time-resolution observation. In addition to quiescent observations, many observations have been 
made selectively towards t he end of outburst as th e star returns to quiescence specifically to study 
DNO and QPO behaviour (jWarner fc Woudtll2006l ) . and several groups of runs taken on consecutive 
nights allow for the study of persistent or evolutionary phenomena across nights. 



Our VW Hyi data are f r om the UCT archive, and have been previous l y analysed inlWoudt fc Warner 



(l2002aT). IWoudt fc Warner! |2002bf) , IWarner. Woudt fc Pretoriusl (120031 1. IWarner fc Woudtl l|20od ) 
Pretorius. Warner k WoudtJ (|2006l ) and IWarner k Pretoriusl (|2008l ) . Details of the individual ob- 
serving runs for outburst and quiescent observations are given in tables [2] and [3] respectively. All 
observations were taken at the Sutherland site of the South African Astronomical Observatory, 
using the 20" , 40" or 74" telescope as indicated. Data acquisition for ob servations on the 20" an d 
40" telescope up to and including run s5248 used the UCT photometer (jNather k Warner Ifl97lh : 
subsequent runs on this telescope and the 74" were made using the University of Cape Town CCD 
(jO'Donoghudll995l ) operating in frame-transfer mode. Most CCD observations were taken in white 
light; observations made with a B filter are marked with an asterisk. 

The data were pre-whitened by removing the mean and first-order trend. VW Hyi lightcurves 
also frequently show a strong orbital hump, which can dominate the periodogram and wavelet 
scalogram. To avoid this, and enable smaller amplitude periodicities to be examined, we removed 
a sinusoid and first harmonic at the orbital hump frequency. 
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4.2. Power Spectrum Estimation 



Torrence &: Compd ()1998l ) 's statistical detection of significant oscillations in the wavelet scalo- 



gram requires estimation of the periodogra m of the n o ise sp ectrum of the signal, a perennial problem 
in noisy astronomical time series analysis. IVaughanl (|2005l ) suggests fitting a linear function to the 
log-periodogram, since his data show a red-noise power spectrum well modeled by power law (i.e. 
linear in log-space). Flickering also gives our data a r ed noise spectru m, but CV flickering has less 
power at very low frequencies than the X-ray data of IVaughanl (|2005l ); we therefore fit a low-order 
polynomial (generally between third and fifth-order) in log-space, using an OLS fitting programme. 

To check that our assumption of a stationary noise structure for each data set was correct, we 
calculated the best-fit noise-model for sections of each run as well as for the entire run; in all cases 
we found that the models for each section were very similar to the model for the complete time 
series. 



4.3. Analysis 

VW Hyi outbursts can be divided into four gro ups: super outburst, lon g normal outburst, 



medium normal outbursts and short normal outbursts (jWoudt fc Warnerll2002al ) . The RASNZ VW 



Hyi data, spanning 1972-2004, have been used to c onstruct average outburs t profiles (figure [2 



onto which each outburst observation can be placed (jWoudt Warne r 2002a). Each observation is 
thus assigned a phase (column 'Start' in tables [2] and [3]) , relative to the arbitrary T =0 d (column 
'T=0' in tables [2] and [3|) of the template. We have assigned quiescent observations a phase relative 
to T =0 d of the next outburst, as we have several quiescent observations occurring a few days 
before an outburst for which we also have data. 

We computed the Morlet ridged scalogram for each run to see if previously detected QPOs were 
detected by the scalogram. We found that all previously detected QPOs were visible in the Morlet 
ridged scalogram, with power above the 95% confidence level and a single continuous instantaneous 
frequency line of stable or slowly changing frequency. We counted the number of extrema (given 
by the maximal points) which occurred during the time span of the confidence contour, and found 
that while 3 cycles (6 maximal points) were sometimes sufficient to have been identified as a QPO 
in the periodogram (e.g. run s0030), usually 5 or more cycles were present. We thus defined a QPO 
as at least 5 cycles (10 maximal points in the MHAT) at power above the 95% confidence level 
(on the Morlet spectrum), or 6 significant cycles spread over two confidence contours at similar 
frequency. We only deviated from this definition when there was a clear triplet of QPOs with 
fundamental, first and second harmonic, and one of these contained less than 5 cycles, in which 
case it was identified as a QPO (e.g. run sl307). We used the instantaneous frequency to find the 
(average) period of the detected QPO. 
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Fig. 27. — Averaged outburst profiles for the four different types of outburst seen in VW Hyi. 
Solid line indicates superoutburst, dashed line indicates long normal outburst, dash-dot indicates 
medium normal and dotted indicates short normal outburst. 
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Table 2. VW Hyi: Outburst observations from the UCT archive. Runs marked with an asterisk 

were taken using a B filter. 



Run Telescope Length Int Outburst HJD Start T = T Start 
(hrs) time (s) Type +2440000 (d) (d) 



s0018 


20 


2.04 


2 


L 


1 572 52 


1572.75 


-0.23 
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1 573 44 


1 579 75 

1UI Li . I U 


69 


s0026 

oUUZiU 


20 


3.32 


5 


L 


1574.50 


1572.75 


1.75 


s0030 


20 
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5 
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1575.47 


1572.75 
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2 


Super 


1662 30 
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1677.20 
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5 
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5 
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-0.83 


s0127 


40 


3.77 


4 


Super 


1677.29 


1677.20 


0.09 


s0129 


40 


1.86 


5 


s 


1691 30 

J- U U -L . OU 


1690 20 

-L U O U . Li U 


1.10 


s0480 


40 


1.73 


5 


kJ LL IJ VI L 


201 7 38 

LiKJ ± I .JO 


2022 65 

Li\J Li Lj .\j\J 


-5.27 


<<0484 

Q U t: O t: 


40 

4:U 


3 91 


4 




2023 31 


9099 65 

LiVJ Li Li .\JO 


66 


<?1 977 

o_L Li 1 1 


40 

4U 


1 78 

1. 1 O 


4 


M 


2352 39 

LiOOL> .Oct 


9354 90 

ZjOUt: . Li U 


-1 81 

l.Ol 


sl307 


40 


1.33 


4 


M 


2354 39 


2354 20 

_ kj \j t: . Zj U 


0.19 


sl322 


40 


4.91 


4 


M 


2355 39 


2354 20 


1.19 


sl571 


40 


1.99 


4 


Super 


2402 32 


2403 40 


-1.08 


sl616 


40 


2.27 


3 


Super 


2404.31 


2403 40 


0.91 


s2241 


40 


2.63 


5 


Super 


2758.41 


2770 00 

Zj 1 1 U . U U 


-11.59 


Szz4o 


a n 
4U 


1.00 


F 




Super 


07cn /in 
z < oy.4U 


z/ /U.UU 


1 n £n 
-1U.DU 


s2623 


40 


3.60 


5 


M 


3515.30 


3514.60 


0.70 


s2911 


40 


3.14 


4 


L 


4934.43 


4937.10 


-2.67 


s2915 


40 


0.73 


4 


L 


4937.28 


4937.10 


0.18 


s2917 


40 


2.18 


4 


L 


4939.49 


4937.10 


2.39 


s3078 


40 


1.95 


4 


Super 


5323.34 


5328.50 


-5.16 


s3416 


40 


1.19 


2 


M 


5967.55 


5967.00 


0.55 


s5248 


40 


4.66 


5 


M 


8202.36 


8201.75 


0.61 


s6138J 


40 


2.09 


4 


M 


11898.28 


11897.85 


0.43 


S6138JI 


40 


3.90 


5 


M 


11898.43 


11897.85 


0.58 


s6184 


40 


1.68 


5 


L 


11957.28 


11957.25 


0.03 


s6316 


40 


5.72 


5 


L 


12354.24 


12353.95 


0.29 


s6528J 


74 


5.81 


4 


L 


12520.22 


12519.65 


0.57 


S6528JI 


74 


4.18 


5 


L 


12520.47 


12519.65 


0.82 


s7222J* 


74 


3.61 


4 


Super 


13007.31 


13006.40 


0.91 


S7222JI* 


74 


3.03 


4 


Super 


13007.46 


13006.40 


1.06 


s7301 


40 


2.94 


5 


M 


13087.24 


13086.10 


1.14 



We give a detailed discussion of several of the runs analysed, showing how DNOs and QPOs 
appear in real data and highlighting the additional informati on which the Mor let ridged scalogram 
makes available. Discussions of further runs are available in iBlackman I ((2008) . 



4.3.1. sOUO, s0122 and s0127 



These three runs are part of a series of runs spanning a superoutburst. Run sOHO was taken 
on during the rise to supermaximum (T =14.90 d), while runs s0122 and s0127 were taken on 
consecutive nights during decline, from T = —0.83 to 0.09 d. 

s0127 (figures [291 and |28|) exhibits a plethora of phenomena, and is an excellent example of both 
the broad frequency range covered by DNOs and QPOs (from periods of 15 s to 1250 s) and the 
intermittent behaviour of these phenomena which have necessitated the development of the Morlet 
ridged scalogram. Figure [281 clearly shows the 27 s DNOs, increasing in period over the duration of 
the run to 34 s. There are also suggestions of the 15 s first harmo nic for the first 3000 s, a nd again 
near the end of the run. Figure [29] shows the 450 s QPO found by I Warner Woudtl (120061 ) . as well 
as evidence of a 1250 s oscillation, changing to 1075 s after 8000 seconds, which we categorize as a 
subharmonic of the 480 s oscillation. 



Runs sOHO and s0122 are included to contrast the Morlet ridged scalograms of data with and 
without significant oscillations. Run sOHO (figure l30"j) shows no oscillations. Although there are 
regions of significant power, the instantaneous frequency fluctuates hugely, and hence these regions 
are not identified as indicating the prese nce of QPOs. In c o ntrast, in s0122 (figure l3"Tj) we see in 
addition to the 410 s QPO identified in IWoudt fc Warner! (|2002al ) a 909 s oscillation persisting 
throughout most of the run, especially prominent in the last half. 



Table 2 — Continued 



Run Telescope Length Int Outburst HJD Start T = T Start 
(hrs) time (s) Type +2440000 (d) (d) 



s7311* 74 5.41 4 M 13139.19 13138.95 0.24 

s7342 40 3.07 5 M 13155.19 13154.90 0.29 

s7621 40 4.03 6 L 13463.23 13463.10 0.13 
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4.3.2. s0129 

s0129 (figure |32|) is a good example of the wavelet spectrum's use in indica ting oscillations of 
l ow coh erence. In addition to the DNO (either DF or lpDNO) at 93 s identified bv lWoudt Warner 
(|2002al ). there is evidence of power at about 33 s throughout the scalogram as well as three runs of 
QPOs with periods 578 s, 658 s and 714 s. The most coherent run lasts for only 5 cycles. Neither 
the 33 s nor the QPOs are visible in the power spectrum. 



4.3.3. sl616 

Run sl616 was taken at the end of return to quiescence after a superoutburst, at 0.912 d, and 
shows a typical triplet of QPO harmonics. A 694 s QPO appears twice in sl616, as well as a 1282 s 
QPO for the first half, and 7 cycles of 333 s QPO, giving the triplet of QPO fundamental, first and 
second harmonic. There is also evidence of an 87 s oscillation near the end of the run (not shown 
in figure |33|) ; this may be either an lpDNO or the DNO fundamental. 



4.3.4. s2241 

Run s2241 (figure [34"|) was taken dur ing superoutburst , and sh ows unclassified QPOs with 



a period of about 700 s, first detected by IWoudt Warnerl (|2002al ). The scalogram reveals the 



previously unrecognized fact that these QPOs are only present during the superhumps - in the top 
panel of figure [3H we show the raw lightcurve of s2241, showing the superhumps, which can be 
matched with the QPOs in the scalogram below. The scalogram enables us to see clearly where 
the oscillations begin and end, and for how many cycles they last. There is some evidence for a 
first harmonic at about 350 s, but this may be due to the non-sinusoidal profile of the QPOs. 

Table |4] summarizes the details of each oscillation detected in our analysis; we include the 
run number, the type of outburst, the start and end times of each oscillation (in seconds from the 
beginning of the run), the period of the oscillation at start and finish (in seconds), and the average 
amplitude (in magnitudes). Oscillations showing non-monotonic period changes were analysed in 
monotonic sections. For oscillations with periods >100 s the number of cycles observed is included. 
DNOs are classified as fundamental (DF), first(Dl) or second harmonic(D2) or lpDNOs (DL). DNOs 
with periods between 70 and 120 s may be either fundamental period DNOs, or lpDNOs; we have 
decided to categorize all DNOs in this period range as lpDNOs. DNO-related QPOs are classified 
as fundamental (QF), first(Ql) or second harmonic (Q2) or subharmonic(QS). QPOs that do not 
appear to be DNO-related are classified as QU, and QPOs from quiescent light curves are classified 
as QQ. The second last column, 'Prev?', indicates whether the oscillation had been detected in 
previous analyses ('y') or not ('n'), and the final column includes comments. 
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Table 3. VW Hyi: Quiescent observations from the UCT archive 



Run 


Telescope 


Length 


Int 


HJD Start 


T = 


T Start 


HJD Next 






(hrs) 


time (s) 


+2440000 


(d) 


(d) 


+2440000 


s0073 


40 


2.83 


5 


1648.25 


1613.00 


-28.95 


1677.20 


s0077 


40 


4.11 


5 


1602.44 


1580.00 


-10.56 


1613.00 


s0085 


40 


1.67 


5 


1604.44 


1580.00 


-8.56 


1613.00 


s0093 


40 


3.33 


5 


1605.45 


1580.00 


-7.54 


1613.00 


s0102 


40 


1.31 


5 


1657.33 


1613.00 


-19.87 


1677.20 


s0105 


40 


1.80 


2 


1660.29 


1613.00 


-16.91 


1677.20 


sl414 


40 


2.89 


4 


2384.30 


2370.00 


-19.10 


2403.40 




Time (s) Amplitude (mag) 



Fig. 28. — Prewhitened lightcurve, Morlet ridged scalogram and periodogram of run s0127. 




Fig. 29. — Prewhitened lightcurve, Morlet ridged scalogram and periodogram of run s0127 for the 
frequency range - 0.01 Hz, showing QPOs. 
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Fig. 30. — Prewhitened lightcurve, Morlet ridged scalogram and periodogram of run sOUO. No 
oscillations are identified. 
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Fig. 31. — Prewhitened lightcurve, Morlet ridged scalogram and periodogram of run s0122, showing 
QPOs. 
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Fig. 32. — Prewhitened lightcurve, Morlet ridged scalogram and periodogram of run s0129. 



- 42 - 




Time (s) Amplitude (mag) 



Fig. 33. — Prewhitened lightcurve, Morlet ridged scalogram and periodogram of run s!616. 



-43- 




1000 2000 3000 4000 5000 6000 7000 8000 9000 0.05 

Time (s) Amplitude (mag) 



Fig. 34. — Raw lightcurve, Morlet ridged scalogram and periodogram of run s2241; note that the 
QPOs visible in the scalogram are only present during superhumps. 
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4.4. Results 



We have identified some 62 new QPOs to add the existing 44, showing that the Morlet ridged 
scalogram is a useful, objective tool for the detection of quasi-periodicities. Figure [35] shows all the 
QPOs detected in our analysis (for ease of viewing, DNO-related QPO harm onic s are shown at the 
i mplie d fundamental). Previously studied QPOs (IWoudt Warnerl (j2002al ) or IWarner fc Woudt 
(|2006l )) are indicated by filled shapes; open shapes indicate QPOs that have been detected using 
wavelet analysis. Circles represent DNO-related QPOs, triangles those QPOs which have no DNO 
counterpart and squares indicate quiescent QPOs. Note that the large numbers of QPOs observed 
around <T(d)< 1.15 is at least partly due to the observing campaign of Brian Warner and Patrick 
Woudt, who selectively observed VW Hyi during late outburst. 



All DNOs previously detected by I Woudt fc Warnerl (|2002al ) and lWarner fc Woudtl (|2006l ) were 
detected in the Morlet scalogram, as well as 2 new DNOs. This is unsurprising: DNOs in general 
have higher coherenc e than QPOs, and a re henc e mo re likely to be detected in periodogram analysis 
such as that used in IWoudt k, Warner! (j2002al ) a nd IWarner fc Woudt] (120061 ) . Out of 24 outburst 
runs showing DNOs, including all those covered in Warner &: Woudtl (|2006 ) . we found DNO-related 
QPOs in each run. In addition, after finding a triplet of QPOs in sl616, close inspection of the FT 
and wavelet showed evidence of DNOs, although of very low coherence. 

In previous work, DNOs with periods > 70 s have been classified as DFs or lpDNOs almost 
arbitrarily; we have chosen to classify every DNO with period > 70 s as a lpDNO. Using this 
classification, we find 8 previously observed lpDNOs, and we have discovered a further 7. 
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Fig. 35. — All QPOs analysed. Circles indicate DNO-related QPOs, triangles indicate unclassified 
outburst QPOs and squares indicate quiescent QPOs. Filled shapes are QPOs previously analysed 
in Warner I, II or IV; open shapes are QPOs added by this study. 



Table 4. Table of QPOs detected in VW Hyi. 



Name 


Type 


Osc 
start 


Osc end 
end 


Cycles 


P start 


P end 


Mag 


Class 


Prev? 


Comments 


s0018 


L 


650 


1900 


4 


346 


346 


0.008 


QF 


n 




s0018 


L 


1000 


2000 




75 


75 


0.007 


DL 


y 




s0018 


L 


1500 


2500 




20 


20 


0.009 


DF 


y 




s0018 


L 


3600 


4600 


3 


286 


286 


0.007 


QF 


n 




s0019 


L 










30 


32 


0.02 


Dl 


y 




s0019 


L 


1600 


3200 


3 


487 


487 


0.05 


Ql 


y 




s0019 


L 


3500 


6250 


7 


370 


400 


0.04 


Ql 


n 




s0019 


L 


5300 


10000 




20 


20 


0.02 


D2 


y 




s0019 


L 


5300 


10000 


4 


510 


510 


0.05 


Ql 


y 




s0026 


L 








12 


1042 


1042 


0.04 


QU 


y 




s0030 


L 


10000 





3 


1124 


1124 


0.04 


QU 


y 




sOlll 


Super 


150 


1300 


4 


307 


307 


0.007 


QU 


n 




s0122 


Super 





5000 


12 


409 


409 


0.005 


QU 


y 




s0122 


Super 


2000 


8000 


7 


909 


909 


0.003 


QU 


n 




s0127 


Super 





2400 




27.9 


28.7 


0.02 


DF 


y 




s0127 


Super 


1000 


7000 


5 


1250 


1250 


0.03 


QS 


n 




s0127 


Super 


2000 


5000 


7 


454 


454 


0.02 


QF 


y 




s0127 


Super 


3600 


4800 




29.3 


29.3 


0.02 


DF 


y 




s0127 


Super 


5100 


8200 




29.7 


30.5 


0.03 


DF 


y 




s0127 


Super 


6000 


8200 


5 


498 


498 


0.02 


QF 


y 




s0127 


Super 


8000 


11000 


3 


1075 


1075 


0.03 


QS 


n 




s0127 


Super 


9100 


12000 




32.2 


32 


0.015 


DF 


y 




s0127 


Super 


10000 


12000 


4 


487 


487 


0.037 


QF 


y 




s0129 


S 





2500 




33.5 


33.5 


0.02 


D2 


n 




s0129 


S 





2500 


4 


578 


578 


0.02 


Ql 


n 


Implied first harmonic 


s0129 


s 


2000 


3500 




93 


93 


0.03 


DL 


y 


DF in Warner - may be DL 


s0129 


s 


3000 


5000 


3 


658 


658 


0.02 


Ql 


n 


Implied first harmonic 


s0129 


s 


3500 


5000 




93 


93 


0.019 


DL 


y 


DF in Warner - may be DL 



Table 4 — Continued 



Name 


Type 


Osc 
start 


Osc end 
end 


Cycles 


P start 


P end 


Mag 


Class 


Prev? 


Comments 


s0129 


S 


5000 


8000 




33 


33 


0.02 


D2 


n 




s0129 


S 


5000 


8000 


4 


714 


714 


0.02 


Ql 


n 


Implied first harmonic 


s0484 


Super 





4000 




30 


32 


0.02 


Dl 


y 




s0484 


Super 





3000 


2.5 


1282 


1282 


0.05 


QF 


y 




s0484 


Super 





12000 


5 


2857 


2757 


0.05 


QS 


n 




s0484 


Super 


2000 


6000 




68 


68 


0.009 


DL 


y 




s0484 


Super 


4000 


6500 




37 


37 


0.02 


Dl 


y 




s0484 


Super 


4900 


7500 


6 


488 


488 


0.03 


Ql 


y 




s0484 


Super 


5000 


14000 


7 


1325 


1325 


0.05 


QF 


y 




s0484 


Super 


6500 


12000 




34 


34 


0.02 


Dl 


y 




s0484 


Super 


8000 


14000 


11 


530 


530 


0.03 


Ql 


y 




sl307 


M 





1000 


2 


500 


500 


0.01 


QF 


n 




sl307 


M 


750 


3000 




31 


31 


0.005 


DF 


n 




sl307 


M 


750 


3000 


5 


250 


250 


0.005 


Ql 


n 




sl307 


M 


1500 


3000 




18 


18 


0.005 


Dl 


y 




sl307 


M 


3500 


4000 




18 


18 


0.005 


Dl 


y 




sl307 


M 


3500 


4500 




38 


38 


0.005 


DF 


n 




sl307 


M 


3500 


4500 


5 


175 


175 


0.01 


Q2 


n 




sl322 


M 





4000 




31 


31 


0.01 


D2 


y 




sl322 


M 





2500 


11 


357 


357 


0.05 


Q2 


n 




sl322 


M 


3000 


6000 


4 


476 


476 


0.05 


Q2 


y 




sl322 


M 


6000 


10000 




29 


30 


0.01 


D2 


y 




sl322 


M 


6000 


10000 


8.5 


465 


465 


0.05 


Q2 


y 




sl322 


M 


10000 


12000 




32 


32 


0.01 


D2 


y 




sl322 


M 


12000 


16000 




28 


31 


0.01 


D2 


y 




sl322 


M 


12500 


15200 


6 


420 


420 


0.02 


Q2 


y 




sl571 


Super 








5.5 


1170 


1170 


0.005 


QU 


y 




sl616 


Super 





2500 


5 


694 


694 


0.03 


Ql 


n 


Implied fundamental 



Table 4 — Continued 



Name 


Type 


Osc 


Osc end 


Cycles 


P start 


P end 


Mag 


Class 


Prev? 


Comments 






start 


end 
















sl616 


Super 





5000 


4 


1282 


1282 


0.05 


QF 


n 


Implied fundamental 


sl616 


Super 


6000 


8000 




87 


87 


0.01 


DL 


n 





Table 4 — Continued 



Name 


Type 


Osc 
start 


Osc end 
end 


Cycles 


P start 


P end 


Mag 


Class 


Prev? 


Comments 


sl616 


Super 


6000 


8000 


7 


333 


333 


0.02 


Q2 


n 


Implied fundamental - low co- 

1 T~\l\T/~\ ""11" " 

nerence DiNOs visible m peri- 
odogram 


s2241 


Super 





2000 


3 


714 


714 


0.02 


QU 


y 




s2241 


Super 


5500 


9400 


5 


740 


740 


0.02 


QU 


y 




s2243 


Super 










75 


75 


0.002 


DL 


y 




s2623 


M 





2000 




22 


22 


0.015 


D2 


y 




s2623 


M 





2000 


5 


408 


408 


0.02 


Ql 


y 




s2623 


M 


1000 


2500 


6 


244 


244 


0.03 


Q2 


y 




s2623 


M 


2000 


4000 




25 


25 


0.014 


D2 


y 




s2623 


M 


4000 


6000 




24 


24 


0.014 


D2 


y 




s2623 


M 


4000 


9000 


3 


1000 


833 


0.005 


QF 


y 


Implied fundamental 


s2623 


M 


5000 


7000 


7 


263 


263 


0.036 


Q2 


y 




s2623 


M 


7000 


10000 


13 


253 


253 


0.035 


Q2 


y 





Table 4 — Continued 



Name 


Type 


Osc 
start 


Osc end 
end 


Cycles 


P start 


P end 


Mag 


Class 


Prev? Comments 


s2623 


M 


7000 


10000 


6 


500 


500 


0.02 


Ql 


n 


s2623 


M 


7500 


10000 




26 


27 


0.021 


D2 


y 


s2623 


M 


10000 


12500 




25 


27 


0.019 


D2 


y 


s2623 


M 


10000 


12500 


4 


416 


416 


0.007 


Ql 


n 


s2915 


L 


700 


1500 


4 


250 


250 


0.006 


Ql 


n 


s2915 


L 


1500 


2500 




20 


20 


0.011 


Dl 


y 


s2915 


L 


1800 


2300 


2 


250 


250 


0.006 


Ql 


n 


s3078 


Super 





7000 




14 


14 


0.004 


DF 


y 


s3078 


Super 


400 


6000 




90 


90 


0.002 


DL 


n 


s3078 


Super 


2000 


4000 


4 


183 


183 


0.004 


QF 


n 


s3078 


Super 


5500 


6500 


4 


182 


182 


0.002 


QF 


n 


s3416 


M 





1000 




25 


25 


0.017 


Dl 


y 


s3416 


M 





1000 


4 


270 


270 


0.026 


Q2 


n 


s3416 


M 


1200 


2500 




27 


27 


0.015 


Dl 


y 


s3416 


M 


1200 


2600 


4 


350 


350 


0.06 


Ql 


y 


s3416 


M 


2600 


4200 




27 


27 


0.017 


Dl 


y 


s5248 


M 





4000 




36 


36 


0.02 


Dl 


y 


s5248 


M 










96 


96 


0.02 


DL 


n 


s5248 


M 


500 


4000 


6 


454 


454 


0.08 


Ql 


n 


s5248 


M 


4000 


8000 




40 


40 


0.005 


Dl 


y 


s5248 


M 


4000 


8000 


5 


694 


694 


0.07 


QF 


n 


s5248 


M 


5500 


13000 


5 


1030 


1030 


0.05 


QF 


y 


s5248 


M 


5500 


13000 


3.5 


2000 


2000 


0.06 


QS 


y 


s5248 


M 


8000 


16000 




40 


40 


0.005 


Dl 


y 


S6138-I 


M 





3000 




25 


25 


0.01 


Dl 


y 


S6138-I 


M 





2200 




25 


26 


0.01 


Dl 


y 


S6138-I 


M 





2200 


7 


317 


317 


0.02 


Ql 


n 


S6138-I 


M 





3000 


3 


1030 


1030 


0.03 


QS 


n 



Table 4 — Continued 



Name 


Type 


Osc 
start 


Osc end 
end 


Cycles 


P start 


P end 


Mag 


Class 


Prev? 


Comments 


s6138-I 


M 


3200 


5000 


3 


667 


667 


0.05 


QF 


n 




s6138-I 


M 


3300 


5000 




26 


28 


0.01 


Dl 


y 




s6138-I 


M 


3300 


5000 


3 


370 


370 


0.06 


Ql 


n 




s6138-I 


M 


5000 


7300 




28 


28 


0.005 


Dl 


y 




s6138-II 


M 





3000 




28 


30 


0.005 


Dl 


y 




s6138-II 


M 





5000 




70 


70 


0.02 


DL 


y 




s6138-II 


M 


5000 


14000 




30.6 


30.6 


0.007 


Dl 


y 




s6138-II 


M 


5000 


14000 


13 


344 


392 


0.04 


Ql 


y 




s6138-II 


M 


6000 


14000 


5 


1450 


1450 


0.02 


QS 


n 




s6138-II 


M 


7000 


14000 


11 


769 


769 


0.04 


QF 


n 




s6184 


L 





1000 




24.8 


27.8 


0.02 


DF 


y 




s6184 


L 





1800 


5 


370 


370 


0.01 


QF 


n 




s6184 


L 


1000 


2300 




26 


26 


0.014 


DF 


y 




s6184 


L 


1800 


2300 




83 


83 


0.01 


DL 


y 




s6184 


L 


2300 


3000 




26 


26 


0.026 


DF 


y 




s6184 


L 


2300 


3800 


3.5 


418 


418 


0.024 


QF 


n 




s6184 


L 


3000 


3500 




25 


25 


0.008 


DF 


y 




s6184 


L 


3500 


3800 




26 


26 


0.014 


DF 


y 




s6184 


L 


3800 


5000 




27 


27 


0.017 


DF 


y 




s6184 


L 


3800 


5000 




74 


74 


0.016 


DL 


y 




s6184 


L 


5000 


5500 




26 


26 


0.013 


DF 


y 




s6184 


L 


5000 


6000 


2.5 


390 


390 


0.01 


QF 


n 




s6184 


L 


5500 


6000 




27 


27 


0.011 


DF 


y 




s6316 


L 





2500 




20 


20 


0.009 


Dl 


y 




s6316 


L 





2000 


6 


294 


294 


0.001 


Ql 


y 




s6316 


L 


2500 


4000 




20 


20 


0.011 


Dl 


y 




s6316 


L 


4000 


6000 




20.5 


20.5 


0.006 


Dl 


y 




s6316 


L 


4000 


6000 




77 


83 


0.008 


DL 


n 





Table 4 — Continued 



Name 


Type 


Osc 
start 


Osc end 
end 


Cycles 


P start 


P end 


Mag 


Class 


Prev? 


Comments 


s6316 


L 


4000 


5000 


5 


227 


227 


0.015 


Ql 


y 




s6316 


L 


6000 


10000 




21.8 


21.8 


0.01 


Dl 


y 




s6316 


L 


10000 


18000 




21 


23 


0.016 


Dl 


y 




s6316 


L 


10000 


15000 




87 


87 


0.012 


DL 


n 




s6316 


L 


11000 


18000 


21 


375 


375 


0.009 


Ql 


y 




s6316 


L 


15000 


18000 




88 


88 


0.011 


DL 


n 




s6528-I 


L 





6000 




18 


23 


0.015 


D2 


y 




s6528-I 


L 





6000 


13 


666 


666 


0.024 


QF 


n 


Implied fundamental 


s6528-I 


L 


6000 


12000 




21 


21 


0.01 


D2 


y 




s6528-I 


L 


6000 


12000 


5 


1299 


1299 


0.024 


QS 


n 


Implied fundamental 


s6528-I 


L 


12000 


18000 




23 


23 


0.011 


D2 


y 




s6528-I 


L 


12000 


18000 


6 


885 


885 


0.013 


QF 


n 


Implied fundamental 


s6528-I 


L 


18000 


20800 




23 


23 


0.009 


D2 


y 




s6528-I 


L 


18000 


20800 




32 


32 


0.008 


Dl 


y 




s6528-I 


L 


18000 


20800 


2 


1052 


1052 


0.017 


QF 


n 


Implied fundamental 


s6528-II 


L 





6000 




24 


24 


0.007 


D2 


y 




s6528-II 


L 





6000 


21 


370 


370 


0.015 


Q2 


n 




s6528-II 


L 


6000 


10000 




24 


24 


0.008 


D2 


y 




s6528-II 


L 


6000 


10000 




37 


37 


0.005 


Dl 


y 




s7222-I 


Super 





2000 




27 


27 


0.01 


D2 


y 




S7222-I 


Super 





10000 


8 


1282 


1282 


0.05 


QF 


n 


Implied fundamental 


S7222-I 


Super 


2000 


4000 




26 


26 


0.005 


D2 


y 




S7222-I 


Super 


2000 


4000 


5 


625 


476 


0.02 


Ql 


n 




S7222-I 


Super 


3000 


6000 




27 


29 


0.011 


D2 


y 




S7222-I 


Super 


6000 


7500 




31 


31 


0.011 


D2 


y 




S7222-I 


Super 


6000 


8000 




38 


38 


0.013 


Dl 


y 




S7222-I 


Super 


6000 


7500 


6 


270 


270 


0.01 


Q2 


y 




S7222-I 


Super 


7000 


9500 


5 


588 


588 


0.04 


Ql 


n 





Table 4 — Continued 



Name 


Type 


Osc 
start 


Osc end Cycles 
end 


P start 


P end 


Mag 


Class 


Prev? Comments 


s7222-I 


Super 


7500 


9500 


29 


29 


0.015 


D2 


y 


s7222-II 


Super 





2500 


34 


34 


0.013 


D2 


y 


s7222-II 


Super 





2500 3 


549 


549 


0.045 


Ql 


n 


s7222-II 


Super 


2500 


4000 


29 


29 


0.015 


D2 


y 


s7222-II 


Super 


4500 


8000 3 


1136 


1136 


0.038 


QF 


n 


s7222-II 


Super 


5000 


9000 8 


333 


333 


0.035 


Q2 


y 


s7222-II 


Super 


5000 


10000 5 


555 


555 


0.037 


Ql 


n 


s7222-II 


Super 


5500 


7500 


28 


28 


0.022 


D2 


y 


s7222-II 


Super 


7500 


8500 


30 


30 


0.016 


D2 


y 


s7222-II 


Super 


8500 


9000 


32 


32 


0.009 


D2 


y 


s7222-II 


Super 


9500 


1000 


34 


34 


0.01 


D2 


y 


s7301 


M 





2000 


33 


33 


0.009 


D2 


y 


s7301 


M 


2000 


4000 


32 


30 


0.015 


D2 


y 


s7301 


M 


4000 


6000 


32 


34 


0.015 


D2 


y 


s7301 


M 


4000 


7000 8 


454 


454 


0.036 


Q2 


n 


s7301 


M 


6000 


7000 


29 


29 


0.014 


D2 


y 


s7301 


M 


7000 


9000 


32 


32 


0.012 


D2 


y 


s7301 


M 


7000 


9000 5 


456 


456 


0.05 


Q2 


n 


s7301 


M 


9000 





29 


31.8 


0.014 


D2 


y 


s7311 


M 








90 


90 


0.014 


DL 


n 


s7311 


M 





1000 3 


312 


312 


0.014 


Ql 


y 


s7311 


M 


1800 


3500 8 


294 


294 


0.017 


Ql 


y 


s7311 


M 


2000 


3500 


19 


19 


0.008 


Dl 


y 


s7311 


M 


3500 


6000 


20 


20 


0.006 


Dl 


y 


s7311 


M 


6000 


15000 


14 


14 


0.015 


D2 


y 


s7311 


M 


6000 


7400 


20 


20 


0.007 


Dl 


y 


s7311 


M 


6000 


10000 12 


316 


316 


0.008 


Ql 


y 


s7311 


M 


7000 


9000 4 


476 


476 


0.028 


QF 


n Implied fundamental 



Table 4 — Continued 



Name 


Type 


Osc 
start 


Osc end 
end 


Cycles 


P start 


P end 


Mag 


Class 


Prev? Comments 


s7311 


M 


7400 


10000 




20 


21 


0.006 


Dl 


y 


s7311 


M 


10000 


13000 


9 


333 


333 


0.028 


Ql 


y 


s7311 


M 


10500 


12500 




21 


21 


0.01 


Dl 


y 


s7311 


M 


12500 







21 


22 


0.009 


Dl 


y 


s7311 


M 


13000 


16000 


6 


345 


345 


0.015 


Ql 


y 


s7311 


M 


13000 


16000 


4 


666 


666 


0.016 


QF 


y 


s7342 


M 










20 


20 


0.006 


Dl 


y 


s7342 


M 





2000 


7 


236 


236 


0.009 


Q2 


y 


s7342 


M 


1800 


6000 


5 


714 


714 


0.006 


QF 


n 


s7342 


M 


2000 


6000 


10 


321 


321 


0.01 


Ql 


y 


s7342 


M 


6000 







14 


14 


0.015 


D2 


y 


s7342 


M 


6000 


8000 




38 


38 


0.009 


DF 


y 


s7342 


M 


6000 


10000 


14 


136 


136 


0.005 


Q2 


n 


s7342 


M 


8000 


9200 


5 


256 


256 


0.015 


Ql 


v 


s7342 


M 


8000 


10000 


4 


500 


500 


0.009 


QF 


y 


s7621 


L 





600 




28 


28 


0.007 


DF 


y 


s7621 


L 


600 


2000 




27 


30 


0.011 


DF 


y 


s7621 


L 


4000 


6000 


5 


454 


454 


0.012 


QF 


y 


s7621 


L 


6000 


11000 


11 


454 


454 


0.027 


QF 


n 


s7621 


L 


8000 


10000 




31 


33 


0.008 


DF 


y 


s7621 


L 


11000 


14000 




18 


18 


0.01 


Dl 


y 


s7621 


L 


11000 


14000 




38 


38 


0.013 


DF 


y 
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5. Discussion and Conclusion 

Lightcurves that contain multiple non-stationary components cannot be satisfactorily analysed 
using the periodogram alone, as the periodogram obscures the time-varying nature of the signal. 
Time-frequency representations allow signals to be viewed in the time and frequency domains 
simultaneously, enabling non-stationary components to be analysed. Various TFRs have been 
used in astronomy, but only for signals that have very high signal-to- noise ratios, cover a narrow 
frequency band, and have only a few components. For signals such as those found in cataclysmic 
variable stars, which contain multiple quasi-periodicities over a broad frequency range, against 
high background noise, it is not clear which of the TFRs used in the literature (if any) is most 
appropriate. 

We have investigated the time- frequency resolution of six TFRs, using three synthetic signals 
mimicking the key features of our data: broad frequency range, multiple components (some close 
together in frequency), amplitude modulation and intermittency. Interference terms are produced 
in many TFRs, so for our data we had to tune the TFRs to completely attenuate ITs, which 
could otherwise be mistakenly identified as signal terms. We found that the CWD could not be 
tuned to completely remove ITs for signals with components that overlap in time, and hence was 
not appropriate for our data. The Gabor spectrogram cannot provide simultaneously high time 
and frequency resolution over the wide frequency range of our data. While both the CKR and 
the ASPWVD gave good time and frequency resolution over a range of frequencies, they could 
not completely attenuate ITs at very high frequencies; the only TFR with achieved this was the 
wavelet scalogram. However, at high frequencies the wavelet scalogram, when tuned for high time 
resolution, did not have good frequency resolution. By computing the wavelet ridges of the Morlet 
scalogram, however, this shortcoming was addressed. We used the maxim points of the Mexican 
Hat scalogram to enable us to detect the extrema of our data. 

We also reviewed the statistical properties of the six TFRs in the presence of noise, and found 
that the wavelet scalogram was the minimum variance estimator for affine TFRs in the presence of 
non-white noise, with a statistical distribution easily computable from the power spectrum. 

We analysed archival data of VW Hyi with the Morlet ridged scalogram, and compared our 
results with those previously published. We found that the Morlet ridged scalogram could be used 
to reliably detect oscillations of only a few cycles (<5) (eg sl307, s6138 and s6184) and we used this 
feature to identify QPOs objectively against criteria based on the appearance of QPOs identified 
in previous analyses. We added some 62 QPOs to the existing 44 of VW Hyi using this method. 

For low amplitude coherent oscillations at short periods (ie < 50 s), it is no surprise that 
the periodogram is still the best means of detection. However, for oscillations in this range that 
have phase changes or gaps (e.g. s0129, sl307, s6528), the Morlet ridged scalogram can be used 
to determine the structure of the oscillation, but only for amplitudes > 0.005 mag. The Morlet 
ridged scalogram can also be used to detect short period oscillations of low coherence, which are 
not visible in the periodogram (e.g. s0129). 
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Oscillations with periods greater than 800 s (eg s0127, s0484, sl322, s6138) fall in a region of 
the periodogram that often shows considerable noise in CV data, making differentiation between 
signal and noise difficult, especially if the amplitude of the oscillation is small, and they have low 
coherence. However, the statistical properties of the scalogram enabled us to fit a model to the 
noise continuum and hence detect significant periodicities. 

The Morlet ridged scalogram also provided a clear picture of the time evolution of different 
frequency components, providing insights unavailable from the raw lightcurve or periodogram. For 
example, groups of spikes in the periodogram were resolved as either a single oscillation with discrete 
period changes (e.g. s0019), or as multiple simultaneous oscillations of similar period (e.g. s2915, 
s2623, s34f6 and sl322). While the QPOs in s2241 had been observed previously, the wavelet 
scalogram showed that they only occur during superhumps, which was not previously recognised. 
We also observed further links between DNO and QPO behaviour, such as period tracking and 
simultaneous initiation/termination of oscillations (e.g. s2623). 



5.1. Software 



The figures of all TFRs in this paper except wavelet scalograms were generated using the 
Matlab Time-Frequency Toolbox, developed by Frangois Auger, Olivier Lemoine, Paulo Gongalves 
and Patrick Flandrin. It is available at http : / /tftb . nongnu . org / \ and is also compatible with GNU 



Octave (available at http://www.gnu.org/software/octave/). An Ansi C version is available at 



http: / /www-lagis.umv-lillel.fr/^davy/toolbox/Ctftbeng.html| 

For the wavelet scalograms, we used the Matlab toolbox of Torrence fc Compo (jl998l ) as a 
starting point (available at http://paos.colorado.edu/research/wavelets/software.html), but have 
written additional instantaneous frequency algorithms and confidence contour functions. All soft- 
ware is available on request. 
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